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Abstract Sersic (1968) generalized the de Vaucouleurs law which follows the projected (observed) 
one dimensional radial profile of elliptical galaxies closely and Dehnen (1993) proposed an analytical 
formula of the 3-dimensional light distributions whose projected line profile resembles the de Vau- 
couleurs law. This paper is involved to recover the Dehnen model and generalize the model to account 



■ for galaxy elliptical shapes by means of curvilinear coordinate systems and employing a symmetry 

principle. The symmetry principle maps an orthogonal coordinate system to a light distribution pat- 



tern. The coordinate system for elliptical galaxy patterns turns out to be the one which is formed 
by the complex-plane reciprocal transformation Z = 1/W. The resulting spatial (3-dimensional) light 
distribution is spherically symmetric and has infinite gradient at its centre, which is called spherical- 
nucleus solution and is used to model galaxy central area. We can make changes of the coordinate 
' system by cutting out some column areas of its definition domain, the areas containing the galaxy 

' centre. The resulting spatial (3-dimensional) light distributions are axisymmetric or triaxial and have 

^ I zero gradient at the centre, which are called elliptical-shape solutions and are used to model global 

elliptical patterns. The two types of logarithmic light distributions are added together to model full 



elliptical galaxy patterns. The model is a generalization of the Dehnen model. One of the elliptical- 
shape solutions permits realistic numerical calculation and is fitted to all R-band elliptical images 
from the Frei et al. (1996)'s galaxy sample. The fitting is satisfactory. This suggests that elliptical 
O \ galaxy patterns can be represented in terms of a few basic parameters. 

' keywords Methods : Analytical ~ Galaxies : Structure 



1 Introduction 

Compared to disc galaxies, elliptical galaxies are more "clean" objects for people to test models on the 
galactic-scale gravitational systems because they have smoother texture of the photographic images. 
Their light distributions can be approximately the linear representation of their mass distributions. 
Real elliptical galaxies are 3-dimensional distributions of light. The observed 2-dimensional elliptical- 
shaped light distributions (i.e., surface brightness) result from the integration of the light along the 
line-of-sight, i. e., the projection on to the 2-dimensional sky plane. Therefore, a complete understand- 
ing of elliptical galaxy patterns is the analytical formulation of their 3-dimensional light distributions, 
i. e., the analytical deprojection of the observed 2-dimensional profiles. The 3-dimensional light distri- 
butions have two aspects: one is the law of their one dimensional radial density profile and the other 
is their 3-dimensional shapes. 

Fortunately, the projected one dimensional radial profile is studied very clear. It decreases smoothly 
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from the nucleus, closely following the Sersic law (1968) 



S,(i?) = a,exp(c7oi?'/") 



(1) 



where as{> 0), cto(< 0) are constants and R is the projected radial distance from the galaxy centre 
along the semi-major axis (however, the spatial (3-dinicnsional) radial distance from the galaxy centre 
is denoted, from now on, by r). The Sersic law, which fits galaxy inner parts very well, indicates 
that galaxy light gradients at their centres are — oo. The Hubble (1930) law suggests an unrealistic 
zero-gradient at the galaxy centres (flattened light distribution). Much attention has been given to 
the deprojection of Sersic law, that is, the proposition of the spatial (3-dimensional) light distributions 
(or equivalently, in the case of spherically symmetric light distributions, the proposition of the spatial 
one dimensional radial profile) whose projected line profiles follow the Sersic law. However, the 
deprojection of the Sersic law is not tractable analytically (Young 1976). An analytical formula of 
the 3-dimensional light distributions whose projected line profile resembles the i?^/" Sersic law or the 
R^^^ de Vaucouleaurs (1948) law is of great interest (Jaffe 1983; Hernquist 1990; Ciotti 1991; Dehnen 
1993; Trijillo 2002). The Dehnen (1993) model is about spherically symmetric light distributions. 



where a, b are constants. 

Another aspect of elliptical galaxy study is their shapes. Elliptical galaxies gain their name by their 
projected apparent shapes. However, their 3-diniensional shapes, axisymmetric or triaxial, are hard to 
infer (Binggeli 1980; Benacchio & Galletta 1980; Binney & de Vacouleurs 1981; Lambas, Maddox, &; 
Loveday 1992; Tenjes et al. 1993; Ryden 1996). The assumption of triaxial shapes is extensively used. 
Stark (1977) proposed a model of triaxial ellipsoid isophotes of constant axis ratios whose projected 
isophote curves are exact ellipticals of constant axis ratios too. The isophote curves of real images, 
however, have radially varying axis ratios. Carter (1978) 's investigation indicates that ellipticity of 
the isophotes is an increasing or peaked function of radius. Analysis of model galaxies suggests a 
similar result (Madejsky &; MoUenhoff 1990; Ryden 1991). In addition, the isophote curves are not 
exact elliptical. They are sometimes boxy-shaped and sometimes discy-shaped (Binney & Merrifield 
1998). Because the model of 3-dimensional nested isophote ellipsoids with varying axis ratios involves 
an infinite amount of freedom of degrees, its projection has to be divided into finite layers and roughly 
fitted to images. 

One goal of the main astrophysical extragalactic studies is to find the symmetry behind the expo- 
nential law of disc galaxies and the Sersic law of elliptical galaxies. With the motivation, I proposed 
to use orthogonal curvilinear coordinate lines to model galaxy patterns (He 2003). Let p{x,y) denote 
spiral galaxy light distribution on 2-dimensional x-y plane (spiral galaxy disc plane). The orthogonal 
curvilinear coordinate system (A, /x) is 



where [x, y) is the common Cartesian coordinates. We always use logarithmic light distribution, 



p{r) = 



b 1 
{r + a)4-T 



(2) 



X = x{X, ji) 

y = y{\ij) 



(3) 



f{x,y) = \ap{x,y). 



(4) 
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Its gradient is V/ whose Cartesian components are (/^, fy). However, its components associated with 
the local reference frames of the curvilinear coordinate system are denoted by n(A, /i) and v{\, fi) (in 
the positive directions of dX and dfi respectively). They are also the directional derivatives to f{x,y) 
in the directions of the curvilinear coordinate lines. Our main idea is that u{X, ji) depends only on A, 
u{X,ij) = u{X) which is the directional derivative to f{X,iJ,) along the coordinate line /i = constant. 
Similarly, v{X, fx) depends only on /x, v(X,fx) = v{fj,) which is the directional derivative to /(A,/x) 
along the coordinate line A = constant. That is, spiral galaxy light distribution f{x,y){= In p{x,y)) 
is such that its gradient vector field can be decomposed in some curvilinear coordinate systems with 
corresponding components depending on single coordinate variables. This is called the single- variable 
requirement, or the symmetry principle. The symmetry principle maps an orthogonal curvilinear 
coordinate system to a logarithmic light distribution. He (2005a) found the symmetric generalization 
of the simple orthogonal coordinate system which is formed by the complex-plane transformation 
Z = . The corresponding light distribution turns out to be spiral galaxy disc. Its radial density 
profile is exponential which is consistent with the known exponential law of spiral galaxy surface 
brightness. In He (2005b), a "disturbed" form of the coordinate system is used to model galactic 
bars. In the present paper, I generalize the idea to 3-dimensional light distributions and present a 
model for elliptical galaxy patterns. An orthogonal coordinate system is employed which is formed by 
the complex plane transformation Z = 1/W. The resulting spatial (3-dimensional) light distribution 
is spherically symmetric and has infinite gradient at the centre, which is called spherical-nucleus 
solution and is used to model galaxy central areas. We can make changes of the coordinate system 
by cutting out some column areas of its definition domain, the areas containing the galaxy centre. 
The resulting spatial (3-dimensional) light distributions are axisymmetric or triaxial and have zero 
gradient at the centre, which arc called elliptical-shape solutions and are used to model galaxy global 
elliptical structures. The two kinds of logarithmic light distributions are added together to model the 
full elliptical galaxy patterns. The model is the generalization of Dehnen's model (2) whose projected 
profile resembles de Vaucouleurs' law. One of the elliptical-shape solutions permits realistic numerical 
calculation and is fitted to all R-band elliptical images from the Prei et al. (1996)'s galaxy sample. 
The model fits real galaxy shapes satisfactorily. This suggests that elliptical galaxy patterns can be 
represented in terms of a few basic parameters. 

Section 2 is devoted to constructing the spherical-nucleus solution. Section 3 presents axisymmetric 
and triaxial elliptical-shape solutions. The two types of solutions are combined to give a full model of 
elliptical galaxy patterns in section 4. A simple discussion of the ellipticity of the projected isophotes 
is given in the same section. Section 5 is the fitting steps of the model to real galaxy images. Section 
6 is simple discussion. 

2 A Model of Galaxy Central Areas: the Spherical- nucleus Solution 

1. Orthogonal coordinates. I propose a spherically symmetric solution of 3-dimensional light distri- 
butions. Spherically symmetric patterns are also axisymmetric: 

f{r) = f{x,y,z) = f{R,z), (5) 

where R = ^x^ ^y^ and the z-axis is the axis about which the 3-dimensional light distribution is 
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Figure 1: The orthogonal 4-circle coordinate system. The {R,z) plane is any one in 3-dimensional 
space which contains the Cartesian z-axis. Therefore, R = ^J{x^ + y^)- 

rotationally symmetric, i. e. axisymmetric. Therefore, we first introduce the standard description 
of orthogonal curvilinear systems defined on a 2-dimensional plane (i?, z) which is any plane 
containing the z-axis (see Fig. 1). A curvilinear coordinate system defined on the plane labels 
each point (i?, z) with an ordered set of two real numbers (A, ji) over a region S on (A, /x) plane. 



The positive increment dA > defines the positive direction of the coordinate line /i = constant. 
Similarly for the coordinate line A = constant. In terms of curvilinear coordinates X = t^, fi = t'^, 
the clement of distance dl between two adjacent points {R, z) = (t"*^, r^) and {R + dR, z + dz) = 
(r^ + dr^^r'^ + dr^) is given by the quadratic differential form (Korn & Korn, 1968), 



R = RiX,fi), 
z = z{X,n). 



(6) 
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dl^ = dR^ + dz'^ = J2J2 gik{r^y)dT'dT^ 



(7) 



i=l fe=l 




(8) 



Orthogonal coordinate system is such that the array {gik) is diagonal, 



dR dR dz dz 



(9) 



We frequently need to trace quantities along curvilinear coordinate lines especially when we draw 
a pattern. Here we consider arc-lengths measured along curvilinear coordinate lines. We denote 



the arc-length along the line /i = constant by s(A, /i) and along the line A = constant by i(A,/x). 
According to the formula (7), they can be calculated as follows, 



ds = ^/gI{dX = ^R'^ + z'^dX = PdX, 
dt = ^/g^dn = ^R'^ + z'^dfjL = Qdfi 

where we have introduced notations for arc-length derivatives, 

p = s'x = ^f^+^ = V9ru 



(10) 



/ (11) 

Q = t'^ = ^R'^ + z'^ = ^. 



2. Com,plex Plane Reciprocal Transformation and Symmetry Principle. Now we present the coor- 
dinate system on which all the models proposed in the paper are based, 

i? = A/(A2+/.2), 

z = ^x/{x^+^,^), (12) 

< A < -|-00, —DO < fl < +00. 

The orthogonal coordinate system can be produced by the complex-plane transformation Z = 
1/W where Z = R + iz and W = X — i/j,. The corresponding arc-length derivatives P, Q are 

P(A,/.) = 4 = l/(A2+^2)^ 

Qix,^^) = t'^ = l/{x^ + ^i^)^p ^ ^ 

The inverse formulas of the coordinate transformation are 

A = i?/(i?2 + z^), 

iJi = z/{R^ + z% (14) 
< i? < -|-oo, —DO < Z < +00. 

Now we look for the logarithmic galaxy light distribution f{R,z) defined on the {R,z) plane. 
The curl of the gradient vector field of the light distribution f{R, z) must be zero. We work on 
curvilinear coordinates and we know that the components of the gradient associated with the 
coordinate are u and v. The zero-curl condition is given in standard mathematical books, e.g. 
Korn & Korn (1968), 

^^{u{\^x)P)-^{v{\^i)Q)=Q (15) 

where P = ^/gii, Q = ^Jg22- Because we require tt, v be functions of single variables A and 
\i respectively (symmetry principle), the above general zero-curl condition changes into (He, 
2005a), 

u{X)P^ - v{fi)Qx = 0. (16) 

Now we need to apply the equation to our circle-shaped coordinate system. Substitution of the 
formulas (13) into above equation, we have the gradient vector field {u,v) for our circle-shaped 
coordinate, 

u{X) = hX,v{ii) = hjjL (17) 
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where h is constant. The galaxy hght distributions along all orthogonal coordinate lines can be 
found by performing path integrations of the following formulas 

df = uds = u{X)PdX, ^^g^ 
df = vdt = v{fi)QdiJ, 

along the coordinate lines /x = constant and A = constant respectively. Choosing the galaxy 
origin to be i? = 0, z = and performing the path integration, we can have our light distribution 
on the {R, z) plane, / = (/i/2)ln(A^ + /U^). We can prove the formula by verifying the rela- 
tion between the physical components of its gradient, u,v, and the contravariant components, 
df /dX,df /dfi. The relation is, u = (df/dX)/P,v = {df /d^)/Q. Using the inverse formulas of 
the coordinate transformation (14), we finally have our logarithmic and direct light distributions 
respectively, 

f{R, z) = {h/2) ln(A2 + ^2) = (/,/2) in 

)V2, (19) 

< i? < +00, —GO < Z < +00. 

where po is another constant. We choose h> because light p = po exp f ^ when R,z^ +oo. 

A Model of Galaxy Central Areas: the Spherical-nucleus Solution. We rotate the 

above (R, z) plane about the z-axis (i.e. the line i? = 0) to achieve our 3-dimensional galaxy light 
distribution which is rotationally symmetric, i. e. axisymmetric. The resulting light distribution 
is 

fn{x, y, z) = {h/2) In 

p„(x,y,z) =po(^.+^)"/' =Po^ (20) 
— oo < x,y,z < +00. 

Wc sec that it is not only rotationally symmetric but also spherically symmetric. The light flux 
at the galaxy centre is +oo and the gradient of the light flux at the same centre is — oo. The 
solution is called spherical-nucleus solution and is used to model galaxy central areas. In fact, 
it corresponds to the first factor of the Dehnen model (2) with 7 being replaced by h. We can 
project the 3-dimensional distribution on any plane to see its surface brightness. We choose the 
x-y plane to be the sky plane. The surface brightness is 

'^{x,y) = J^^ p{x,y,z)dz 

= '^Po!n^^f"dz 

= 2poS^'-{^f'^dz (21) 

= m) 

_ r((fe-i)/2) 

where Gamma function is used and the integral can not be carried out to be an elementary func- 
tion of the parameter h. The surface brightness at the galaxy centre is +00 which is inconsistent 
with the Sersic law. 



3 Models of Galaxy Global Patterns: the Elliptical-shape Solutions 

In this section I present axisymmetric and triaxial models of elliptical galaxy patterns. 
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Figure 2: The orthogonal coordinate system is obtained by simply eutting out two infinite eolumns 
of areas in Fig. 1. The two columns are parallel to and centreing on the Cartesian i?-axis and z-axis 
respectively. The widths of the two columns are 2i?o and 2zo respectively. 



1. A model of 3-dimensional axisymmetric light distributions. The model is a corrected 
form of the above 3-dimensional light distribution (20), 

fisix, y, z) = {h/2) In (^^,^^)4(^^^|^|), , 

pis{x, y, z) = /^o( (fi„+fi)4(.o+k|)^ )''^' ^^^^ 

< i? < +00, —DO < Z < +00 

where Rq{> 0) and zq{> 0) are constants. This is equivalent to cutting out an infinite long 
cylinder of light distribution (of radius i?o) around the z-axis and cutting out an infinite plane 
layer of light distribution (of width 2zq) parallel to and centreing on x-y plane from the spherically 
symmetric distribution (20). The corresponding orthogonal coordinate system is demonstrated 
in Fig. 2 which is the one in Fig. 1 with corresponding parts being cut out. The corresponding 
coordinate equation system is 

R = A/(A2 + ^2) _ 

\z\ = \^i\/{\^ + ^,'')-zo, (23) 

< A < i?o(-Ro + 4)hh < < (i?o + 4)1 

We can see that the 3-dimensional light distribution (22) is rotationally symmetric and its 
projection on any sky plane is apparently elliptical except on the x-y plane. 

Now, we project the 3-dimensional light distribution on the sky plane which is a plane containing 
the X-axis and is obtained by rotating the x-y plane about the x-axis with an angle (the 
inclination angle of the axisymmetric galaxy, see Fig. 3). The coordinates {x' , y') on the sky plane 



y' 



Figure 3: we project the 3-dimcnsional light distribution on the sky plane which is a plane containing 
the X-axis and is obtained by rotating the x-y plane about the x-axis with an angle 60 (the inclination 
angle of the axisymmetric light distribution) . 

(see Fig. 3) and the coordinate z' along the line-of-sight are related to the original coordinates 
{x,y,z) as follows 

y = y' sec 9q + z' sin ^o, (24) 
z = z' cos ^0- 

The surface brightness on the sky plane after the projection is given as the following integration, 

nx\y')=!^^p{x,y,z)dz' 

= On f+°°( 1 )'^/'^dz' "^^^^ 

The x'-axis is the major or minor axis of the 2-dimensional elliptical light distributions S(x', y'). 
The integral can not be carried out to be an analytic formula. It can be numerically calculated. 
However, it can be analytically solved in some special cases. 

Firstly, take h = 3,9q = 0, and zq = 0. The corresponding surface brightness is 

^{x',y')=2po/iRo + RT (26) 

where R' = ^ x'"^ + y'"^. This is exactly the Hubble (1930) law. Secondly, take /i = 3, = 0) and 
zo 7^ 0. The corresponding surface brightness is 

= (^^^^ - 7(?T(^TW)^ ^''^ 

which is a generalized Hubble law. The general case {h > 3) requires numerical calculation. 
However, the meaning of the parameters is simple. The model of axisymmetric light distributions 
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(22) involves four parameters pQ,h,RQ,ZQ. The central light density of the light distribution 
is proportional to po while h represents galaxy light concentration and all other parameters 
{Rq, zq and xq, yo iii the following part (v)) affect the elliptical shapes of the 3-dimensional light 
distributions. The projection of the light distribution on to the sky plane (the integration (25)) 
can not be expressed by the F function of h as we did in (21). The total luminosity, i.e. the 
summation of the light distribution through whole space is finite when /i > 3. 

2. The 3-parameter axisymmetric model. Our model of axisymmetric light distributions 
(22) involves 4 parameters {po^h^RQ, zq). Its projection on to the sky plane (the integration 
(25)) can not be expressed by the V function of h as we did in (21). However, if we choose Rq 
to be zero, our 4-parameter model changes into a 3-parameter model (poi h,zo), 

P2s{x, y, z) = po{^—^-^———)f'/^ (28) 
+ + [zo + \z\y 

and its projection on to the sky plane can be evaluated by the F function of h and the integrals 
on finite intervals, 

= Po(/o^°°( ..^,.4„^.). f/^d/ 

= Poim^Tik^)'^' + iwHtT^)''')dz' (29) 

where 

R± = V^'^ + (y'T^osin0o)' ^3Q^ 
Z± = ±y' tan + zq cos 

are not dependent on the integration variable z' . Finally, our 3 parameter model gives light 
distributions on the sky plane as follows, 

E(x',2/')=Po( 

W^'^^ r{h/2) + Jz+/r+{t+^) (31) 

"^1 /V;F r((fe-l)/2) rO ( 1 \h/2j„\\ 

+ ^f^\-^ IW2) ^ iZ-/R-yT+^) "SJJ- 

This is exactly the model which is used in section 5 to be fitted to real elliptical galaxy images. 

3. Orthogonal coordinates and zero-curl equations in the case of 3-dimensional space. The curvi- 
linear coordinate system in 3-dimensional space is 

X = x{X, fx, u), 

y = y(A,p,z^), (32) 

Z = z{X, jJL, v). 

The arc-length derivatives are, 

P = s\ = ^{x'^ + y'^ + z'^) = ^^, 

Q = t'= Vi^u + y'' + z'^) = yfm. (33) 



„/2 I „,/2 I J2\ 



G = V(< + y[/ + 4') = v^. 
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We work on curvilinear coordinate system and the components of the gradient vector of the 
logarithmic light distribution /(x, y, z) associated with the coordinate system are u, v and w. 
Because they are functions of the single variables A, /x, v respectively, the zero-curl condition 
consists of the following equations, 

u{X)P^ - v{fi)Qx = 0, 

u{X)P, - w{u)Gx = 0, (34) 
v{lJi)Qv - w{v)G^ = 0. 

Now we present the 6-sphere coordinate (A, /x, u) which is the generalization of the plane 4-circle 
coordinate (formulas (12) and Fig. 1) to the case of 3-dimensional space, 

= A/(A2 + ^2 ^ 1.2)^ 
y = /x/(A2 + + 1/2)^ 

Z = i./(A2 + ;,2 + ^2)^ ^ ) 

— GO < X, fJ,,!^ < +00. 

In Section 2, the plane 4-circle coordinate lines are rotated about the 2;-axis. The resulting 3- 
dimensional coordinate surfaces consist of one family of 2-spheres, one family of toroids, and the 
other family of all planes containing the z-axis. Here the coordinate system has triaxial symmetry 
and the coordinate surfaces consist of three families of 2-spheres. Therefore, the coordinate is 
called a 6-sphere coordinate (Moon & Spencer, 1961). The corresponding arc-length derivatives 
P, Q, G are 

P(A,/X,Z.) =l/(A2+^2 + ^2) ^^^^ 

= Q(A,^,z^) = G{X,n,v). 
The inverse formulas of the coordinate transformation are 

A = x/(x2 + y2 + 2;2), 
/X = y/(x2 + y2 + z2), 

V = z/(x^ + y + z ), 
—GO < x,y,z < +00. 

The zero-curl equation system (34) determines the gradient vector field {u, v, w) of our 3- 
dimensional logarithmic light distribution, 

n(A) = hX, v{n) = h/j,, w{v) = hv (38) 

where /i is a constant. The galaxy light distributions along all orthogonal coordinate lines can 
be found by performing path integrations of the gradient components, 

fn = (h/2) ln(A2 + ^2 ^ ^2) ^ In 

Pn{x,y,z) = po(^^+^)'^/2 = p^i. (39) 
—GO < x,y,z < +00. 

We choose h > because light p = poexpf when x,y,z +oo. The light distribution 
turns out to be the spherical-nucleus solution. However, we can achieve triaxial and other light 
distributions based on the result, using the cut-out method. 



10 



4. A model of spherically symmetric light distribution with zero gradient at the galaxy 
centre. We can cut out a sphere (radius a) of light distribution centred at the galaxy nucleus. 
The resulting light distribution is 

f3s{x,y,z) = {h/2) 

Pzs{x,y,z) = poj^:^, (40) 
< r < +00 

In fact, it corresponds to the second factor of the Dehnen model (2) with 4 — 7 being replaced 
by h. 

5. A model of 3-dimensional triaxial light distributions. We can cut out three infinite plane 
layers of light distributions parallel to and centreing on the three Cartesian coordinate planes 
respectively from (39). The three planes are the x-y plane, y-z plane, and z-x plane. The widths 
of the three layers are 2xo,2yo,2zo respectively. The resulting logarithmic light distribution is 

f4s{x, y, z) = {h/2) In (^g+|^|)a+(2^o+|2/l)'+(2o+kl)^ ' 

P4s{x, y, z) = Po((^Q+|^|)^+(J,o+|J/|)^+(^o+|^|)0''/^ (41) 
—00 < x,y,z < +00 

where xo{> 0),yo(> 0) and -2;o(> 0) are constants. 

Now, we project the 3-dimensional light distribution on the sky plane to see the surface bright- 
ness. Take the x-y plane to be the sky plane as an example. To achieve an analytic result, take 
^ = 3. The resulting surface brightness is 

(I , 2^ ^ 

^ Vi^'^+{^o+\x\r+{yo+\y\r)' 

which is similar to the formula of the circularly symmetric light distribution (27). But here the 
surface brightness is non-circularly symmetric even when the sky plane is the x-y plane. 



A Model of Full Galaxy Patterns and Some Symmetry Results. 

1. A Model of Full Galaxy Patterns. To achieve a full account of the central areas as well 
as the global shapes, we need to add the logarithmic light distribution of the spherical-nucleus 
solution in section 2 to one of the logarithmic light distributions of elliptical-shape solutions in 
section 3. The summation of two logarithmic light distributions is equivalent to the multiplication 
of their direct light distributions, PnPis- If) instead of multiplication, we add together the direct 
light distributions, p„ + pis, then we would sec that elliptical galaxies consist of two apparent 
components. Real galaxy images, however, do not show such apparent components. The Dehnen 
model is recovered 

p{r) = Pn{r)p3s{r) = po-^^-^——^ (43) 

if we choose /ii = 7, /t2 = 4 — 7. This gives a model of spherically symmetric patterns. The 
projected radial profile of the model light distribution resembles the de Vaucouleurs law. Its 
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detailed discussion is given in Dehnen (1993). All models PnPis, i = 15 2,3,4, in the present 
paper have similar properties. Here, we pay special attention to the PnP2s axisymmetric model 

p{x,y,z) = PnP2s ,,,, 



which gives account of both galaxy central area and elliptical shape. The model is called full 
3-parameter axisymmetric model. The projection of the light distribution on to the sky plane 
can not be evaluated by the F function of h and the integrals on finite intervals as we did in (31) 
for the original 3-parameter model. 

2. EUipticity as a function of radius. Firstly, let us study the ellipticity of the projected 
isophote curves of the original 3-parameter model p2s- The 3-dimensional light distribution of 
the model is the result of a spherically symmetric one with a layer of light centred on the x-y 
plane being cut out. Therefore, if the sky is a plane containing the z-axis then the projected light 
distribution (surface brightness) on the sky plane has approximately elliptical isophote curves 
(discy-shaped) whose ellipticity is approximately 



The ellipticity decreases monotonically with radius R from e = 1 at the galaxy centre to e = 
at i? = -l-oo. Actually, we see 1 > ei > e{R) > €2 > because p2s has a fiat light distribution 
near the galaxy centre and the image is always taken within a finite area. 

Secondly, we study the full 3-parameter axisymmetric model, p{x,y,z) = pnP2s- The ellipticity 
e(i?) is a peaked function of R changing from e = at the galaxy centre to e = at i? = +oo, 
because the spherical-nucleus solution has dominant logarithmic light density (~ +00) near the 
galaxy centre while all elliptical-shape solutions have finite values (flat light distributions) near 
the galaxy centre. Actually, we see an increasing ellipticity function if zq is large enough, or 
a decreasing function if zq is small enough, or a peaked or flat function if zq has intermediate 
values. These are consistent with Carter (1978) 's image analysis. The ellipticity function and 
the isophote-curve shapes (discy or boxy) change with the orientation of the line of sight and 
with the models of elliptical-shape solutions. 

3. Light-strength invariance. Note that there is always an arbitrary factor po accompanying 
our model light distribution (see, e.g. (21) and (22)). This is because our models involve the 
logarithmic light distributions f{x,y,z) and the symmetry principle deals with their gradients 
only and leaves an integration constant for f{x,y,z). This is called light-strength invariance 
which indicates that bright galaxies can share the same patterns with dim galaxies. Because 
galaxy visible mass distributions are approximately proportional to their light distributions, 
p{x, y, z) can be mass distribution or light distribution in our model. In our fitting program, 
however, the projection of p{x,y,z) on to the sky plane is fitted to the light distributions on 
galaxy images. 

The satisfactory fitting of our model to real images (see section 5) suggests that elliptical galaxy 
patterns can be represented in terms of a few basic parameters given in the model. Real galaxies 




1 1 



(44) 



(a;2+j,2_|_(_j,g + |^|)2)h2/2 




(45) 
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are further constrained by their physical process. For example, the shape of the radial profile and 
the shape of isophote curves change with luminosity of the objects. Our mathematical model 
gives larger freedom of degrees and bright galaxies are allowed to share the same patterns with 
dim galaxies. 

4. Scale invariance and the unit of lengths. We can multiply a constant A to the coordinate 
equations (12) to form a new coordinate system (scale transformation), 

Then we have P = AP,Q = AQ,u = u/A, and v = v/A. The light distribution determined by 
the new coordinate system is 

p{R,z) = p{R,z). (47) 

This is called scale invariance. The scale transformation is equivalent to choosing different units 
for the Cartesian coordinates R, z. The scale invariance says that our method of generating 
galaxy patterns does not depend on the choice of length units. 

5 Fitting the 3-Parameter Axisymmetric Model to Elliptical Galaxy 
Images 

Only the projections of the 3-parameter axisymmetric model p2s a^nd the triaxial model p^s are found 
to be expressions of T functions and integrals on finite intervals. The projections of other models 
need further investigation. The direct numerical calculation of the integrals on infinite intervals takes 
significant computer time to obtain a reasonable result. 

In the present paper we fit the 3 parameter axisymmetric model (31) to real elliptical galaxy 
images and we find that the fitting is satisfactory except the areas near the galaxy centres. To achieve 
satisfactory fitting of full range of images, we need to apply the full model introduced in section 4 
whose projection, however, is not an expression of T functions and integrals on finite intervals and 
needs deep investigation. 

Fitting Step 1: Rotating the real elliptical galaxy image. We need to fit the expression S(a;', y') (31) 
to real elliptical galaxy images. Because we cut out an infinite layer of light distribution centreing on 
the x-y plane and we rotate the light distribution about the x-axis, the x'-axis must be the major axis 
of the projected 2-dimensional elliptical light distributions S(x',y') (see (24) and Fig. 3). Therefore, 
we need to find the galaxy centre of the image and rotate the image about the centre so that the major 
axis of the elliptical galaxy image is in the horizontal direction as the ones of the model patterns are. 

The sky level needs to be subtracted. The resulting light distribution of the galaxy image is 
denoted by 

S(i,j) (48) 

where {i,j) are the image pixel coordinates. 

Fitting Step 2: One more fitting parameter I: the side length of the square variance domain of 
{x',y') of the model formula. The distance d to us of a galaxy can be obtained by the Hubble law 
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Figure 4: The Prei et al.{1996) R-band images, NGC 3377, 3379, 4125, minus our theoretical fitting 
patterns respectively result in the residual images. Explanation of the spots on the fitting images are 
given in the Fitting Step 3. 




Figure 5: The Frei et aL(1996) R-band images, NGC 4365, 4374, 4406, minus our theoretical fitting 
patterns respectively result in the residual images. Explanation of the spots on the fitting images are 
given in the Fitting Step 3. 
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d = cz/H, given the Hubble constant H (e. g., H=75 km/(s-Mpc)) and the redshift z of the galaxy. 
Therefore, the physical position (£, y) on an elliptical galaxy image is known because the angles 
subtended by the positions against us are measurable by the observer. Therefore, real galaxy light 
density T;{x,y) can be considered the function of real physical position {x,y). The analytic formula 
(31), however, is a specific mathematical function E of independent variables x' ,y' which can be 
denoted by whatever letters and have no physical meaning. However, if real galaxy light distributions 
of all galaxies can be expressed by a single universal function then we can say that x' ,y' are 
proportional to physical position variables if we can testify that the formula (31) approaches to be the 
function. The test consists of two steps. The real image is obtained from a finite square range of x, y 
on a galaxy and the side length of the range, I, is called the physical size of the data. The formula (31) 
is defined on — oo < y' < +oo and we need to choose a finite square range of x', y' for the theoretical 
formula, that is, we need to choose a part of the theoretical light distribution (31) which is compared 
to the digital image data. If the theoretical formula approaches to be the exact expression of all galaxy 
light distributions then there must exist the part of the theoretical light distribution (31) which best 
fits the real galaxy data, along with the fitted values of all other parameters the formula involves. The 
side length I of the fitted square range of x',y' is called the fitted size of the data. We calculate the 
ratio of the fitted size to the physical size of the data, ifl. We do this for a sample of galaxies. This 
is the first step. The second step is very simple. We see if all the ratios are approximately identical. 
If they are, we can say that the galaxy light distributions I](x,y) of all galaxies can be expressed by 
a single universal function and the formula (31) approaches to be the function and the mathematical 
variables x',y' have physical position meaning. A physical unit is qualified for the variables x',y', 
denoted by Ch. Its equivalent value to the meter or pc is given by the identical ratio. Testing the 
identical ratio is equivalent to testing the proportionality of galaxy redshifts z to the fitted lengths Iq 
of one arcsecond distances on the galaxy images. One arcsec distance on a galaxy image means that 
the corresponding physical distance on the galaxy makes an angle of one second at us (the observer). 
The fitted size I divided by the total arcseconds of one side of the corresponding galaxy image is the 
fitted arcsecond length Iq. Whenever we are given a mathematical formula we should not lose the 
chance to test if it approaches to be the universal physical one. 

The first step of the test is explained on more detail. As explained in the previous paragraph, the 
length of the square range of x' , y' needs to be considered a fitting parameter (wc do not worry about 
the centre position of the range because the centre of the model light distribution is exactly the galaxy 
centre). We try different values of the length I. For each value, the square domain is divided into 
M X M small squares whose length is l/AI. The value of the theoretical light at the i x j-th square, 
S(i, j), is compared with the image light value at the i x j-th pixel, X!(i, j). The total relative error is 

M 

err= ^ \t{i,j) -T.{i,j)\/t{i,j). (49) 

Note that we take the relative error not the absolute error because both the pixel light values and 
the signal-to-noise ratios are very large near galaxy centres. Even the best fitting function, the Sersic 
law, leaves large absolute errors near the galaxy centres. To reduce the effects of the galaxy central 
areas, I choose relative error for the fitting. The set of fitted values of the above scale I and all other 
parameters of the model formula makes the minimum err. 
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Table 1: The Fitted Parameter Values for the S ersic Law and the 3-Parameter Model 

NGC redshift Iq Gg n gq 



[/ arcsec] 


3377 


0.0022 


0.040 


72782 


3.2 


-5.41 


3379 


0.0030 


0.087 


58621 


3.0 


-3.88 


4125 


0.0045 


0.014 


56853 


3.5 


-6.75 


4365 


0.0042 


0.038 


63469 


3.2 


-4.65 


4374 


0.0035 


0.050 


51938 


3.0 


-4.34 


4406 


-0.0008 


0.017 


47698 


3.3 


-5.58 


4472 


0.0033 


0.036 


67274 


3.0 


-4.34 


4621 


0.0014 


0.037 


26526 


2.7 


-4.35 


4636 


0.0031 


0.017 


50065 


3.6 


-5.58 


5322 


0.0059 


0.033 


134785 


4.0 


-6.42 


5813 


0.0066 


0.036 


34088 


3.0 


-4.34 


NGC 


^0 


Po 


h 


e° 






[Ch/arcsec] 








[Ch] 


3377 


0.15 


462862 


4.90 


79 


5.29 


3379 


0.13 


615071 


5.00 


37 


4.78 


4125 


0.28 


479334 


4.06 


83 


9.51 


4365 


n 2Q 


1 46908 


3 46 


74 


7 33 


4374 


0.35 


159427 


3.50 


38 


6.52 


4406 


0.09 


14994 


3.62 


66 


2.99 


4472 


0.32 


122068 


3.24 


66 


5.71 


4621 


0.12 


241165 


4.86 


71 


4.41 


4636 


0.30 


103418 


3.27 


66 


7.91 


5322 


0.53 


990340 


3.83 


79 


13.21 


5813 


0.13 


547757 


4.75 


65 


5.05 



Fitting Step 3: Run computer program to find the fitted values of all parameters. I wrote a 
computer program for the fitting, that is, I try all possible values of the set of parameters and calculate 
the corresponding (31) and the error (49). Finally I find the fitted values of the parameters 

which make the minimum err (see the result of the second half of the Table 1). The corresponding 
fitted elliptical pattern is shown in the middle columns of Figures 4 and 5) . Note the two symmetric 
small spots on the central vertical direction of each fitted pattern in the Figures. This is due to almost 
zero values of and i?_ at or near the two points (see (30)). Therefore, the integrals in (31) are on 
almost infinite intervals and their numerical calculations result in large errors by the computer which 
runs over a limited time. Now we know that it is important that wc have an analytic model whose 
calculation can be taken on computers in a reasonable amount of time. 

The second half of the Table presents the fitting result for all 11 R-band elliptical galaxy images of 
Frei et al. (1996) 's sample (the sample consists of nearby bright galaxies and the images are cleaned by 
the original authors). The second column of the half of the Table is the fitted lengths Iq of one arcsec 
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0.002 0.004 0.006 0.008 0.01 

redshift z 

Figure 6: The fitted lengths (fiUed cireles) of one-aresee distances on elhptieal galaxy images are 
approximately proportional to the galaxy redshists. This is not true for the fitted lengths (open 
circles) based on one dimensional Sersic law. 

distances on the images. Fig. 6 indicates that the fitted lengths Zq are approximately proportional to 
galaxy redshifts. This is equivalent to say that the ratios of the fitted sizes to the physical sizes of the 
image data are approximately constant. That is, the galaxy light distributions y) of all galaxies 
can be expressed by a single universal function and the formula (31) approaches to be the function 
and the mathematical variables x', y' have physical position meaning. A physical unit is qualified for 
the variables x',y', denoted by Ch. Its equivalent value to the meter or pc is given by the identical 
ratio as follows 

l[Ch] = 0.00027[Mpc] = 270[pc]. (50) 

Figures 4 and 5 demonstrate the 2-dimensional fitting result. The fits along the galaxy major axes 
are shown in Fig. 7. The observed profiles are indicated by the symbols. The solid lines are the fitting 
result taken from Figures 4 and 5. The dotted lines are the fitting result of the Sersic law (1). Similar 
to 2-dimensional model fitting, when fitting a one dimensional Sersic law I need to select a variance 
domain of R for the function (1). However, I consider the length I of the domain {0,1) as a fitting 
parameter. I try different values of the length I. For each value, the domain is divided into M small 
intervals whose length is l/M. The value of the theoretical Sersic light at the i-th interval, T>s{i), is 
compared with the image light value at the i-th pixel along the major axis, S(i). The total relative 
error is 

M 

err = 5^|S(z)-S,(i)|/S(z). (51) 

1=1 

The set of fitted values of the length / and all other parameters of the Sersic law makes the mini- 
mum err. In fig. 7, the value at the upper right corner of each panel is the fitted value of I for the 
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Figure 7: The observed profiles are indicated by the symbols. The solid lines are the fitting result 
taken from Figures 4 and 5. The dotted lines are the fitting result of the Sersic law (1). 

corresponding galaxies. The fitted size I is divided by the total arcseconds of the major axis on the 
corresponding galaxy image. The result is Iq which is the fitted length of one arcsec distance on the 
image. The fitted values of the parameters in the Sersic law together with Zq are presented in the 
first half of Table 1. The values of Iq fitted by the Sersic law do not suggest a linear correlation 
with redshifts (open circles in Fig. 6) . One of the reasons is that an image is the projection of the 
corresponding 3-dimensional galaxy light distribution and the one dimensional Sersic law cannot give 
account of the galaxy shape. 

6 Conclusion 

One goal of the main astrophysical extragalactic studies is to find the symmetry behind the exponential 
law of disc galaxies and the Sersic law of elliptical galaxies. With the motivation, I proposed to 
use orthogonal curvilinear coordinate systems and a symmetry principle to model galaxy patterns 
(He 2003). The symmetry principle maps an orthogonal coordinate system to a light distribution 
pattern. The coordinate system for elliptical galaxy patterns turns out to be the one which is formed 
by the complex-plane reciprocal transformation Z = 1/W. The resulting spatial (3-dimensional) 
light distribution is spherically symmetric and has infinite gradient at its centre, which is used to 
model galaxy central area. By cutting out some column areas of its definition domain the coordinate 
system can also give axisymmetric or triaxial spatial light distributions. The two types of logarithmic 
light distributions are added together to model full elliptical galaxy patterns. The Dehnen model is 
recovered and generalized. The isophote curves of the projected patterns can be discy or boxy and 
their ellipticity changes with radius. One of the elliptical-shape solutions permits realistic numerical 
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calculations and is fitted to all R-band elliptical images from the Frei et al. (1996)'s galaxy sample. 
The fitting is satisfactory. This suggests that elliptical galaxy patterns can be represented in terms of 
a few basic parameters. 

The coordinate systems describing 2-dimensional disc galaxy patterns are based on the coordinate 
system which is formed by the complex-plane exponential transformation (He 2005b). He (2005b) 
found the evidence that the coordinate system can be considered a "free fall" coordinate space and a 
non-geometrized yet relativistic stellar dynamics was developed based on the principle of coordinate 
cancellation of gravity. A simple explanation of constant galactic rotation curves is given in He (2005b) . 
The present paper indicates that elliptical galaxy patterns can be explained by the coordinate system 
(35). Whether the coordinate system cancels gravity in 3-dimensional space for elliptical galaxies 
needs deep investigation and is left to be future work. 
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